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ABSTRACT 


Optimal  maintenance  policies  are  designed  to  reduce  the  number  of  system  failures 
and  minimize  the  cost  of  repair  by  scheduling  planned  replacements.  In  this  area  the 
problem  of  updating  the  maintenance  policy  using  the  past  maintenance  history  has  not 
been  adequately  solved.  In  this  thesis  we  study  a  sequential  estimation  procedure  in  a 
nonparametric  setting  to  estimate  the  age  replacement  policy  that  minimizes  long  run 
expected  maintenance  costs. 

This  thesis  begins  with  the  discussion  of  the  concepts  of  preventive  maintenance,  age 
replacement  policies,  the  settings  of  our  simulation  model,  and  a  detailed  description  of 
the  sequential  estimation  procedure.  We  include  examples  using  actual  replacement  data 
which  demonstrate  the  usefulness  of  the  sequential  estimation  procedure.  Monte-Carlo 
methods  are  used  to  study  the  behavior  of  estimated  optimal  age  replacement  policy  for 
different  sample  sizes,  costs  and  underlying  system  life  distributions.  We  also  make 
comparison  with  Frees  and  Ruppert's  (1985)  sequential  procedure  for  estimating  optimal 
age  replacement  policies.  These  comparisons  show  that  our  sequential  estimation  pro¬ 
cedure  is  competitive  and  for  large  sample  sizes  performs  better  than  the  Frees  and 
Ruppert's  procedure.  Finally,  we  will  introduce  a  graphical  method  to  estimate  the  op¬ 
timal  age  replacement  policy. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Optimal  maintenance  policies  are  designed  to  reduce  the  number  of  system  failures 
and  minimize  the  cost  of  repair  by  scheduling  times  for  replacements.  By  far,  most  of 
the  research  in  this  area  has  been  from  the  modeling  standpoint.  Even  in  the  most  basic 
scenario,  where  the  underlying  system  lifetimes  are  assumed  to  be  independent  and 
identically  distributed,  the  problem  of  updating  the  maintenance  policy  using  the  past 
maintenance  history  has  not  been  adequately  solved.  In  this  thesis  Monte-Carlo  meth¬ 
ods  are  used  to  study  a  particular  nonparametric  procedure  which  updates  estimates  of 
an  optimal  age  replacement  policy  after  each  replacement.  In  particular,  we  will  focus 
on  finding  conditions  under  which  this  sequential  procedure  does  well. 


B.  MAINTENANCE  POLICIES 

The  effectiveness  of  a  working  system  depends  not  only  on  the  innate  properties 
built  into  it  in  the  design  and  production  stages,  but  also  on  the  quality  of  its  operation, 
maintenance  and  repair.  Maintenance  policies  are  designed  to  reduce  the  incidence  of 
system  failure  by  scheduling  maintenance  actions.  For  example,  a  maintenance  policy 
is  warranted  when  the  failure  of  a  system  during  actual  operation  is  costly  or  dangerous 
(i.e.,  when  the  cost  C,  of  unscheduled  maintenance  due  to  system  failure  is  more  than 
the  cost  C2  of  scheduled  maintenance  before  system  failure).  If  the  system  is  character¬ 
ized  by  a  failure  rate  that  increases  with  age,  it  may  be  wise  to  replace  it  before  it  has 


aged  too  greatly.  This  is  the  concept  of  preventive  maintenance.  Broadly  speaking, 


preventive  maintenance  is  the  total  of  all  service  functions  aimed  at  maintaining  and 


improving  reliability  performance  characteristics  and  concerns  itself  with  such  activities 
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as  the  replacement  and  repair  of  systems,  inspections,  testing  and  checking  of  working 
parts  during  their  operation. 

In  the  preventive  maintenance  model  developed  here,  it  will  be  assumed  that  the 
maintenance  action  returns  the  equipment  to  the  “as  good  as  new”  condition,  thus 
providing  the  same  services  as  equipment  that  has  been  “replaced”.  This  assumption 
implies  that  the  times  between  failure  are  independent  and  have  the  same  distributions. 
If  this  is  not  the  case,  the  model  needs  to  be  modified.  In  addition,  the  model  assumes 
that  the  scheduled  and  unscheduled  replacement  costs  C,,  C2  remain  constant  over  time. 

The  two  most  common  replacement  policies  are  the  policy  based  on  age  (age  re¬ 
placement)  and  the  policy  based  on  time  (block  replacement).  Under  age  replacement 
a  sv  stem  is  always  replaced  at  the  time  of  failure  or  at  age  T,  whichever  occurs  first. 
Under  a  block  replacement  policy  the  system  is  replaced  upon  failure  and  at  times  T  , 
2T  ,  3T  ,  ....  ,  etc.  By  its  nature,  age  replacement  is  administratively  more  difficult  to 
implement,  since  the  age  of  the  unit  must  be  recorded.  On  the  other  hand,  block  re¬ 
placement,  although  simpler  to  administer  since  the  age  of  the  unit  need  not  be  recorded, 
leads  to  more  frequent  replacement  of  relatively  new  items  [Ref.  1:  pp.  178-182].  Thus, 
in  this  thesis  we  only  consider  age  replacement  policies. 

We  define  an  optimal  age  replacement  policy  as  the  age  replacement  policy  which 
yields  the  smallest  long  run  expected  replacement  costs  per  unit  time.  To  determine  an 
optimal  age  replacement  policy  we  require  explicit  knowledge  of  the  system's  life  dis¬ 
tribution.  But  in  a  real  world  scenario,  we  may  not  know  the  life  distribution  explicitly, 
thus  the  optimal  age  replacement  policy  needs  to  be  estimated.  Estimation  when  a  fixed 


sample  of  independent  and  identically  distributed  (/./.</.)  system  life  times  were  available 


has  been  examined  in  detail  [Refs.  2, 3, 4, 5],  To  gather  such  i.i.d.  data,  experimental  sys¬ 


tems  must  be  left  in  service  until  failure.  Therefore,  the  experimenter  can  not  implement 
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an  age  replacement  policy  to  achieve  cost  savings  while  collecting  data  for  estimation. 
When  the  luxury  of  observing  system  lifetimes  until  failure  is  not  available,  a  more  cost 
effective  approach  is  to  estimate  sequential^.  By  this  we  mean  that  the  estimator  of  the 
optimal  age  of  replacement  is  updated  after  each  system  replacement.  In  addition,  a 
system  under  observation  is  subject  to  a  replacement  policy  that  is  close  to  the  best  es¬ 
timated  policy  so  Atr.  The  nonparametric  procedure  of  this  type  which  we  will  study  is 
described  in  detail  in  Chapter  2.  The  difficulty  with  such  a  sequential  procedure  is  that 
implementing  an  age  replacement  policy  while  collecting  data  results  in  censored  obser¬ 
vations.  Thus,  there  is  a  trade-off  between  the  goals  of  controlling  cost  and  gathering 
data  for  estimation.  This  conflict  is  particularly  acute  in  the  nonparametric  setting 
where  information  about  the  right  tail  of  the  distribution  needed  to  estimate  the  optimal 
replacement  age  can  not  be  obtained  under  heavy  censoring. 

This  thesis  will  include  small,  moderate  and  large  sample  studies  by  Monte-Carlo 
simulation  to  determine  the  behavior  of  the  sequential  procedure  under  various  condi¬ 
tions.  In  Chapter  2,  description  of  the  age  replacement  policy  and  the  sequential  esti¬ 
mation  procedure  are  given.  Included  in  this  chapter  is  an  example  of  the  sequential 
estimation  procedure  applied  to  actual  data.  Chapter  3  establishes  the  setting  for  the 
simulation  study,  with  a  quick  overview  of  the  results  of  that  study,  A  detailed  analysis 
of  the  simulation  results  is  given  in  Chapter  4.  In  that  chapter,  we  also  present  results 
comparing  our  sequential  estimation  procedure  with  another  sequential  estimation  pro¬ 
cedure.  In  Chapter  5,  we  describe  a  graphical  method  for  estimating  the  optimal  age 
replacement  policy.  Conclusions  and  recommendations  are  given  in  Chapter  6. 


II.  THE  SEQUENTIAL  ESTIMATION  PROCEDURE 


A.  THE  OPTIMAL  REPLACEMENT  AGE 

Let  {X„)  be  an  independent  and  identically  distributed  (i.i.d.)  sequence  of  positive 
random  variables  with  distribution  function  F.  The  sequence  A'„  A\  ....  represents  the 
sequence  of  system  lifetimes  that  would  be  observable  if  the  sv stems  \vere  replaced  at 
failure.  It  is  intuithely  clear  that  it  is  not  advantageous  to  use  an  age  replacement  policy 
for  a  system  whose  lifetime  has  a  decreasing  failure  rate  or  constant  failure  rate  since  this 
would  result  in  a  replacement  by  a  new  system  whose  properties  are  worse  than  or  just 
the  same  as  the  original  one.  Let  C,  be  the  cost  of  an  unscheduled  replacement  (at  sys¬ 
tem  failure)  and  C2  be  the  cost  of  a  scheduled  replacement  (before  system  failure),  where 
C,  >  C2.  Under  an  age  replacement  policy  with  scheduled  replacement  at  age  t,  standard 
results  from  renewal  theory  [Ref.  6:  p.  87J  can  be  used  to  show  that  the  long  run  ex¬ 
pected  cost  per  unit  time  is  given  by 


m  = 


C)  x  F(t)  +  C2  x  5(0 

[  S(.v)  dx 
Jo 


t 


where  S(.t)  =  1  —  F(x)  is  the  survival  function.  The  numerator  is  the  expected  cost  of 
one  replacement  under  the  age  replacement  policy,  and  the  denominator  is  the  expected 
time  between  replacement.  Under  some  fairly  general  conditions,  there  exists  unique  and 
finite  time  <p*  where  /?(/)  attains  a  global  minimum  [Ref.  7:  pp.  161-168).  For  example, 
a  sufficient  condition  for  the  existence  of  <f>*  is  that  Fhave  failure  rate  ?.(x)  that  strictly 
increases  to  infinity.  We  will  assume  through  out  the  sequel  that  </>*  exists. 
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When  such  a  </>*  exists,  it  is  the  optimal  age  of  replacement,  in  the  sense  that  a 
policy  with  age  replacement  at  (j>*  will  have  minimum  long  run  expected  costs  per  unit 
time.  When  </>*  is  unknown  then  it  needs  to  be  estimated  from  system  lifetime  data. 

B.  THE  SEQUENTIAL  ESTIMATION  PROCEDURE 

Let  {</>*}  be  the  sequence  of  estimators  of  </>*  where  the  estimator  (/>*  is  based  on 
data  from  the  first  it  replacements.  If  our  goal  was  just  to  estimate  </>*,  then  we  would 
construct  </>*  based  on  observing  n  complete  system  lifetimes.  In  this  case  all  replace¬ 
ments  would  be  unscheduled  replacements,  at  system  failure,  with  cost  C,.  We  wish, 
however,  to  control  replacement  costs  by  implementing  age  replacement  policies  while 
collecting  data  to  estimate  </>*.  There  is  clearly  a  trade-off  between  gathering  data  for 
estimation  (i.e.,  observing  system,  as  long  as  possible)  and  controlling  costs  by  imple¬ 
menting  the  age  replacement  policy  with  the  best  estimated  scheduled  replacement  age 
so  far.  Therefore,  we  compromise  and  choose  the  scheduled  replacement  age  for  the 
/fth  system  to  be  £n  where  (/>*,  <  c„  <  oo.  Intuitively,  it  seems  reasonable  to  take 
c„  =  </>*,  +  a„  where  {«„}  is  a  sequence  of stants  that  decrease  to  zero  as  n  -» oo.  How 
to  choose  the  sequence  {«„}  is  an  open  question  that  we  study  in  detail  in  the  Monte- 
Carlo  simulations  that  follow. 

Assuming  that  the  (a„)  sequence  and  c:„  the  replacement  age  for  the  first  system 
under  observation,  are  given,  the  procedure  to  compute  the  estimators  { </>*  }  developed 
in  [Ref.  8]  follow's: 

1.  Determine  /th  system's  scheduled  replacement  time  by 

Zi  =  1  +  ai  ■ 
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2.  At  the  /th  replacement  observe  the  system  lifetime  Xt  or  <?,  whichever  comes  first. 
Let  Z,  =  min(^V„  £,)  and  <),  =  I(X,  <  q,)  where  1(A)  is  indicator  function  of  the  set  A. 
In  other  words,  if  system  is  replaced  before  failure,  then  the  replacement  time 
Z,  =  q,  and  S,  =  1,  otherwise  Z,  =  X and  5,  =  0. 

3.  Sort  Z„  Z2, Z,  into  an  ascending  sequence  Za„  Zl2), . Z<0  with  <5(1„  <5(2), . <5(0 

ordered  according  to  the  ordering  of  Z(1),  Z(2>, Z(ir 

4.  Use  the  data  (Z(1),  <5(1)),  (Z(2),  <5(J)),  (Zw,  S(0)  to  estimate  5  s  1  -  F.  The  estimator 

A 

S  of  S  is  defined  as 


*«-  n 

{ / :  2(I>  <  O 


Note  that  although  5  is  formally  identical  to  the  Kaplan- Meier  estimator 
[Ref.  9:  pp.  34-35]  based  on  randomly  right  censored  data,  it's  properties  are  sub¬ 
stantially  different  because  the  underlying  data  pairs  (Z,,  <5,)  ,  (Z2,  S2) .  ,  (Z„  St) 

are  not  independent. 

5.  Use  the  estimator  in  step  4  to  estimate  R(t)  as 


m  = 


C,  x  F(t )  +  C2  x  S(t) 


I’f  A 

S(X) 


dx 


6.  The  minimum  of  the  estimator  R(t)  does  not  exist.  In  addition  if  <5(0  =  0,  i.e.,  if  the 
largest  replacement  age  observed  so  far  occurred  before  system  failure,  then 

A 

K(r)  -*  0  as  n  -+  oo.  To  take  care  of  these  two  problems,  after  /  replacements,  wc 
estimate  by  $*,  where  0*  attains  the  infimum  of  R(t)  over  o<,t<  Z(0.  In  fact 
0*  is  easy  to  calculate.  It  can  be  shown  [Ref.  8],  that  0*  is  equal  to  the  Z,  among 
1  <j<,  i  such  that 

R(Zj)  =  min  {  R(Z^) ,  R(Z^) . . R(Z^ }. 


The  procedure  is  then  repeated. 
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This  sequential  estimation  scheme  means  that  after  (/  —  1)  replacements  the  exper¬ 
imenter  has  available  a  scheduled  replacement  age  c,.  The  replacement  cost  for  the  /th 
system  is  C,  if  Xt  <  otherwise  the  replacement  cost  is  C.  Thus  with  this  sequential 
estimation  procedure  the  actual  total  replacement  cost  for  the  first  n  systems  that  are 
observed  is 


n 

Cn  =  Yj  (C1  X  31  +  C2  X  (  1  ‘  Sl  )  )> 
/=1 

and  the  total  operating  time  for  the  n  systems  is 


n  n 

tn  =  Yj inin  w  ^ ) =  E  z<  • 

1=1  i=i 

Let  t  >  0,  and  a,  >  t  for  /  =  1,  2, .  Then  in  [Ref.  8]  it  is  shown,  that  as  n  ->  oo,  </>* 

converges  to  </>*  with  probability  1.0,  and  the  actual  cost  per  unit  time  ~~  converges 
with  probability  1.0  to  R(<£*).  Since  s  >  0,  can  be  arbitrarily  small,  for  practical  pur¬ 
poses  we  can  choose  a  sequence  { an }  such  that  an  ->  0.  It  is  clear  that  choice  of  {«„}  can 
have  considerable  impact  on  how  well  0*  estimates  0*  and  on  how  large  the  actual 

c 

costs  ~  are.  An  important  feature  of  this  thesis  is  to  use  simulation  to  provide  practical 
guidelines  for  choice  of  {a„}. 

C.  EXAMPLE:  REPLACEMENT  COST  ANALYSIS  OF  THE  TRACTOR  AT 
ENGINE  FAILURE 

To  show  how  the  sequential  procedure  works,  we  consider  the  data 
[Ref.  10:  p.  560],  summarized  in  Table  1,  which  includes  the  age  at  engine  failure  of  22 
brand  new  D9G-66A  Caterpillar  tractors  and  the  calendar  date  of  the  failure. 
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Assume  the  unscheduled  replacement  cost  C,  =  S200.00  and  scheduled  replacement 
cost  C,  =  S  100.00.  If  engines  are  replaced  onlv  at  failure,  with  no  ase  reolacement 

*  w  1  ^  7  V  * 

A 

policy  (i.e.,  c,  =  co  for  /  =  1, 2, .... ,  22),  then  at  each  replace:.  Z,  =  X„  and  F  becomes 
the  empirical  distribution  based  on  the  system  lifetimes  observed  so  far.  Table  2  on 
page  9  gives  the  estimates  of  </>*  calculated  after  each  replacement. 


Table  1.  THE  AGE  OF  TRACTOR  FAILURE 


Sequence  of  failure 

Tractor's  number 

Date  of  failure 

Age  in  hours  when 
failure  occurred 

1 

16 

4-15-65 

5161 

2 

12 

6-2S-65 

5278 

3 

13 

8-1-66 

6378 

4 

15 

8-3-66 

6578 

5 

14 

9-14-66 

6385 

6 

17 

10-26-66 

6717 

7 

19 

4-3-67 

5556 

8 

10 

5-8-67 

2690 

9 

18 

11-1-67 

6869 

10 

11 

3-26-68 

6259 

11 

8 

8-8-69 

4394 

12 

22 

10-31-69 

6150 

13 

2 

4-16-70 

5085 

14 

5 

6-1-70 

6052 

15 

7 

8-10-70 

7774 

16 

9 

9-21-70 

10517 

17 

6 

6-7-71 

6367 

18 

1 

6-16-71 

8230 

19 

20 

7-30-71 

3286 

20 

3 

10-11-71 

3826 

21 

21 

1-21-72 

4815 

22 

4 

5-8-72 

10950 
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From  Table  2,  we  have 


total  replacement  cost  =  22  x  200.00  =  4400.00 , 


total  operating  hours  = 


135299,  and 


....  4400.00  dollars 

cost  per  unit  time  =  •  ,  =  0.032o2  -  — -  . 


Table  2.  ESTIMATION  WITH  NO  AGE  REPLACEMENT  POLICY 


i 

?/ 

A'; 

8, 

z, 

n 

CO 

5161 

1 

5161 

5161 

2 

CO 

5278 

mm 

5278 

5161 

3 

CO 

6378 

mm 

6378 

5161 

mm 

oo 

6578 

MM 

6578 

5161 

5 

oo 

6385 

i 

6385 

5161 

6 

oo 

6717 

mm 

6717 

5161 

7 

oo 

5556 

i 

5556 

5161 

8 

oo 

2690 

i 

2690 

5161 

9 

oo 

6869 

i 

6869 

5161 

10 

oo 

6259 

n 

6259 

5161 

11 

oo 

4394 

■9 

4394 

5161 

12 

oo 

6150 

l 

6150 

5161 

13 

oo 

5085 

wm 

5085 

5085 

14 

oo 

6052 

wm 

6052 

5085 

15 

oo 

7774 

mm 

7774 

5085 

16 

oo 

10517 

wm 

10517 

5085 

17 

CO 

6367 

6367 

5085 

18 

CO 

8230 

i 

8230 

5085 

19 

oo 

3268 

i 

3268 

5085 

20 

oo 

3826 

i 

3826 

5085 

21 

oo 

4815 

in 

4815 

5085 

22 

oo 

10950 

l 

10950 

5085 
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We  now  use  the  sequential  estimation  procedure  described  in  the  last  section  on  the 
same  set  of  data,  with  the  same  costs  C,  =  200.00,  C2  =  100.00.  For  this  example,  let 
c,  =  oo,  so  that  first  replacement  will  occur  at  system  failure  and  let  a,  =  - hours 
for  i  =  2,  3, ....,  22.  The  sequence  {<?„}  used  in  this  example  was  chosen  for  convenience. 
No  special  attempt  was  made  to  find  the  { an }  sequence  that  gave  the  best  results.  Ap¬ 
plying  the  sequential  estimation  procedure  introduced  in  last  section,  the  sequential  es¬ 
timates  of  0*  and  the  data  are  summarized  in  Table  3. 


Table  3.  ESTIMATION  WITH  AGE  REPLACEMENT  POLICY 


i 

«; 

c, 

a; 

s, 

Z, 

mm 

0 

oo 

5161 

1 

5161 

5161 

2 

500.00 

5661.00 

5278 

mm 

5278.00 

5161 

3 

250.00 

5411.00 

6378 

0 

5411.00 

5161 

mm 

166.67 

5327.67 

6578 

0 

5327.67 

5161 

5 

125.00 

5286.00 

6385 

0 

52S6.00 

5161 

6 

100.00 

5261.00 

6717 

0 

5261.00 

5161 

7 

83.33 

5244.33 

5556 

0 

5244.33 

5161 

8 

71.43 

5232.43 

2690 

l 

2690.00 

5161 

9 

62.45 

5223.50 

6869 

0 

5223.50 

5161 

10 

55.56 

5216.56 

6259 

0 

5216.56 

5161 

11 

50.00 

5211.00 

4394 

l 

4394.00 

5161 

12 

45.45 

5206.45 

6150 

0 

5206.45 

5161 

13 

41.67 

5202.67 

5085 

i 

5085.00 

5085 

14 

38.46 

5123.46 

6052 

0 

5123.46 

5085 

15 

35.71 

5120.71 

7774 

0 

5120.71 

5085 

16 

33.33 

5118.33 

10517 

0 

5118.33 

5085 

17 

31.25 

5116.25 

6367 

0 

5116.25 

5085 

18 

29.41 

5114.41 

8230 

0 

5114.41 

5085 

19 

27.78 

5112.78 

3268 

mm 

3268.00 

5085 

20 

26.32 

5111.32 

3826 

i 

3826.00 

5085 

21 

5110.00 

4815 

l 

4815.00 

5085 

22 

23.81 

5108.81 

10950 

0 

5108.81 

5085 
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This  procedure  yields 


total  replacement  cost  = 

total  operating  hours  = 
cost  per  unit  time  = 


E 

1=1 

22 

Z 

/=» 


(200.00  x  <5/  +  100.00  x  (1 -<?,)) 
Z,  =  107395.48,  and 


=  0.02793  ■df11  — 
10/.  9.3.48  hour 


3000.00 , 


At  each  updating,  both  procedures  give  identical  estimates  4>*,  i  -  1,2, 22.  But 

the  final  cost  per  unit  time  when  no  age  replacement  policy  was  implemented  is  larger. 

By  varying  unscheduled  maintenance  cost,  it  can  be  seen  in  Table  4  that  as  the  ratio  of 
C 

—  increases,  the  benefit  of  applying  an  age  replacement  policy  while  gathering  data 
V-2 

from  estimation  becomes  greater. 


Table  4.  COMPARISON  OF  MAINTENANCE  COSTS 


Sched¬ 
uled  re¬ 
placement 
cost 

System 
opera¬ 
tion  time 

Total 
mainte¬ 
nance  cost 

Estimated 
optimal 
replace¬ 
ment  time 

Average 
cost  with 
Age  Re¬ 
placement 
Policy 

Average 
cost  with 
Failure 
Replace¬ 
ment  Pol¬ 
icy 

200 

100 

107395.48 

3000 

5085 

0.02793 

0.03252 

300 

100 

106358.44 

3800 

4815 

0.03573 

0.04878 

400 

100 

102600.13 

4300 

4394 

0.04191 

0.06504 

500 

100 

99185.13 

4600 

3826 

0.04638 

0.08130 

600 

100 

98029.13 

5200 

2690 

0.05305 

0.09756 

700 

100 

96919.44 

5200 

2690 

0.05365 

0.11382 

800 

100 

91040.44 

5700 

2690 

0.06261 

0.13008 

900 

100 

86865.44 

6200 

2690 

0.07137 

0.14634 

1000 

100 

95194.75 

670C 

2690 

0.07038 

0.16260 
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III.  SIMULATION  SETTING 


A.  UNDERLYING  LIFE  DISTRIBUTION 

Guidelines  For  choosing  a  sequence  {a„}  for  estimation  will  depend  on  three  factors, 
the  underlying  distribution  F,  the  costs  C„  C2  and  the  sample  size  A'  where  we  intend  to 
stop  sampling.  The  two  factors  that  we  assume  are  known,  are  the  sample  size  N  and 
the  costs  C,  and  C2.  In  much  of  the  simulation,  we  concentrate  on  fixed  costs  C,  =  5.0 
and  C2  =  1.0.  Other  costs  are  also  considered,  but  in  much  less  detail.  We  choose 
sample  sizes  A’  -  10,  50  and  250  to  reflect  small,  moderate  and  large  sample  sizes.  The 
third  factor  is  the  underlying  system  life  distribution  which  in  general  will  be  unknown. 
For  simulation,  we  choose  this  distribution  to  be  Weibull  with  shape  parameter  a  and 
scale  parameter  ).,  where  the  density  is  given  by 

J[i)  =  «  /.  (/U)a-1  e~  W  for  t  >  0. 

The  Weibull  distribution  has  failure  rate 

;.(*)  - « (/it)*-1 . 

When  a  >  1.0,  the  failure  rate  is  strictly  increasing  to  infinity.  Thus,  for  Weibull  dis¬ 
tributions,  with  a  >  1.0,  a  unique  and  finite  optimal  replacement  age  </>*  exists.  The  ten 
different  Weibull  distributions  used  in  the  simulation  have  a  values  1.1,  1.2, .....  1.9,  2.0. 
This  selection  of  a  values  gives  us  a  range  of  distributions  which  become  more  like  the 
exponential  distribution  as  a  decreases  from  2.0  to  1.1.  To  make  fair  comparisons  be¬ 
tween  Weibull  distributions,  the  scale  parameter  /.  was  chosen  so  that  the  expected  sys¬ 
tem  lifetime  E( X)  =  2.0.  See  Figures  1  gnd  2,  for  plots  of  the  Weibull  densities  and 
corresponding  failure  rates. 
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WEIBULL  DENSITY  FUNCTION 


Figure  1.  The  Weibull  density  function  f(0  with  E(X,)  =  2.0 


i 

I 


Figure  2.  The  failure  rate  of  Weibull  distribution  with  £(X,)  =  2.0 


B.  OPTIMAL  REPLACEMENT  TIME 

When  the  distribution  of  system  lifetime  is  Weibull,  the  expected  long  run  average 
cost  per  unit  time  under  a  simple  age  replacement  policy  with  scheduled  replacement  at 
age  t,  is  given  by 


m= 


c.  ( i  -  <r(;-r)‘)  +  c2*-('1-0, 

1‘e-^dx 

Jn 


See  Figure  3  for  a  plot  of  R(t)  when  the  underlying  life  distribution  is  Weibull  wun  shape 
parameter  a  varying  from  1.1  to  2.0.  For  eacii  curve  on  Figure  3  the  optimal  replace¬ 
ment  time  (p*  can  be  located  on  the  x-axis  at  the  minimum  point. 

Table  5  on  page  16  gives  the  optimal  replacement  times  <p*  for  different  values  of 
the  shape  parameter  a  when  the  unscheduled  replacement  cost  C,  =  5.0,  and  the 
scheduled  replacement  cost  C2  =  1.0.  The  optimal  replacement  times  <f>*  in  Table  5  are 
computed  from  simulation  results.l  Eacii  <p*  was  generated  using  60,000  pseudo  ran¬ 
dom  Weibull  lifetimes.  Included  in  Table  5,  is  the  probability  that  a  system  will  be  re¬ 
placed  before  failure  under  the  optimal  age  replacement  policy, 


P{Xi<(}>*)  =  >\ 

From  Table  5  on  page  16,  we  observe  that  the  optimal  replacement  time  <f>*  in¬ 
creases  as  the  shape  parameter  a  decreases  from  2.0.  For  values  of  o.  close  to  1.0,  the 
performance  of  a  new  system  will  not  differ  greatly  from  the  old  one  which  is  still  in  use. 
In  this  case,  very  little  is  gained  by  replacing  the  system  before  failure  at  the  higher  cost 
C,.  The  larger  values  of  <£*,  insure  that  a  small  percentage  of  replacements  will  be  made 

1  Since  the  optimal  replacement  time  <f>*  comes  from  a  simulation  result,  it  varies  slightly 
with  the  number  of  pseudo  random  variables  used  and  the  seed  numbers  used  to  generate  them. 
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COMPARISON  OF  LONG  RUN  EXPECTED  AVERAGE  COSTS 
(  C,-5.0  ,  C2-1.0  ) 


Figure  3.  Long  run  expected  average  cost  curves  with  E(X()  =  2.0 

before  failure,  which  is  what  we  desire  if  the  system's  life  distribution  is  close  to  expo¬ 
nential. 

C.  OVERVIEW  OF  SIMULATION  RESULTS 

A  wide  range  of  {crn}  sequences  have  b<\»"  selected  for  study.  Their  functional  forms 
are  given  in  Table  6  on  page  18.  Recall  that  q,  (i.e.,  the  scheduled  replacement  time  or 
scheduled  censoring  time  for  the  /th  system  lifetime)  is 


fi- 


+  alf 


/  =  2,  3, . 
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Table  5.  ESTIMATED  OPTIMAL  REPLACEMENT  TIMES  OF  WEIBULL 


MODEL  WITH  E(X.)  =  2.0 


Shape  pa¬ 
rameter  a 

Scale  pa¬ 
rameter  X 

Optimal  replace¬ 
ment  time 

Long  run  expected 
optimal  replace¬ 
ment  cost  R{<{>*) 

P(,Y,  <  <£*) 

2.0 

0.443113 

1.1792297 

1.8024244 

0.238937 

1.9 

0.443682 

1.1846266 

1.8709278 

0.255175 

1.8 

0.444643 

1.2031918 

1.9470497 

0.277008 

1.7 

0.446122 

1.2661085 

2.0291557 

0.315235 

1.6 

0.44S2S7 

1.3114500 

2.1227818 

0.347835 

1.5 

0.451373 

1.4643221 

2.2216024 

0.415706 

1.4 

0.455712 

1.6853428 

2.3185654 

0.49S967 

1.3 

0.461788 

2.3491287 

2.4112816 

0.670973 

1.2 

0.470328 

3.7472544 

2.4730082 

0.861093 

1.1 

0.482456 

12.9815130 

2.4966393 

0.999460 

We  have  arbitrarily  selected  c,  =  1.0  and  let  a,  =  0.0  for  all  {«„}  sequences.  For  each 
a  ,  sample  size  X  and  sequence  {a„}  ,  we  have  simulated  the  results  of  sequentially  esti¬ 
mating  </>*.  Each  simulation  is  based  on  generating  1000  sequences  of  system  lifetimes. 

Let  RjX,  j  =  1,2, .  1000  be  the  actual  replacement  cost  per  unit  time  for  the  first 

N  replacements  of  they'th  repetition  of  a  simulation.  For  each  simulation,  the  perform¬ 
ance  of  the  sequential  estimation  procedure  is  first  evaluated  by  computing 


1000  /n  n/x^\\2 


M  l 


1000 
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where  MSE  is  the  average  squared  difference  of  the  actual  replacement  cost  per  unit  time 
from  the  minimum  long  run  expected  replacement  cost  per  unit  time.  Because  the  ulti¬ 
mate  goal  of  estimating  0*  is  to  reduce  costs,  we  have  chosen  to  evaluate  performance 
by  comparing  at  actual  costs  per  unit  time  to  /?($*),  rather  than  looking  at  the  mean 
square  error  of  the  estimators  of  </>*.  Tables  7,  8  and  9  summarise  the  simulation  results 
of  the  {«,,}  sequences  that  performed  best  for  different  values  of  a  and  sample  sizes  .Y, 
with  fixed  costs  C,  =  5.0,  C2  =  1.0.  The  criteria  for  selecting  the  best  performing  {«„} 
sequence  is  to  choose  the  sequence  which  gives  the  smallest  MSE.  In  cases  where  {«„} 
sequences  had  nearly  equal  MSE,  we  select  the  sequence  with  the  lowest  actual  cost  per 
unit  time  averaged  over  the  1000  repetitions. 


1000 


J=  i 


We  have  also  plotted  the  best  performing  {a.}  sequences  using  the  three  sample  sizes 
of  10,  50  and  250  with  various  values  of  the  shape  parameter  a  in  Figure  4  on  page  22. 
From  these  plots  and  the  results  from  Tables  7,  8  and  9,  we  conclude  that  in  general 
when  a  decreases  from  2.0  (i.e.,  the  underlying  life  distribution  is  becomes  more  expo¬ 
nential),  the  best  performing  {«„}  sequence  tends  to  have  larger  values.  In  the  next 
chapter  we  provide  a  detailed  analysis  of  these  simulation  results,  that  reveals  how  the 
best  performing  { an }  sequences  are  effected  by  sample  size,  underlying  distribution  and 
costs. 


» 
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Table  7.  BEST  SIMULATIOiN  RESULTS  FOR  SMALL  SAMPLE  SIZE  (N=  10) 


Shape 
param¬ 
eter  « 

{«„}  functional  form 

Actual 

average 

cost 

Rs 

MSE  of 
actual 
average 
cost 

Esti¬ 

mated 

optimal 

replace¬ 

ment 

time  <j>* 

MSE  of 

i>* 

2.0 

1 

1.725  +  0.06  x(i-  1) 

2.01557 

0.36262 

1.26406 

0.27861 

1.9 

1 

1.725  +  0.06  x(t-  1) 

2.06978 

0.38469 

1.25512 

0.28013 

1.8 

1 

1.725  +  0.01  x(/-  1) 

2.13750 

0.41103 

1.27617 

0.32854 

1.7 

1 

1.400  +  0.03  x(<-  1) 

2.21857 

0.44410 

1.33542 

0.38985 

1.6 

1 

1.400  +  0.03  x(/-  1) 

2.29646 

0.49618 

1.36779 

0.41863 

1.5 

1 

1 .400  + 0.03  x  (f  -  1) 

2.37696 

0.54815 

1.39615 

0.45544 

1.4 

1 

1.400  +  0.03  >:(/-  1) 

2.47533 

0.64276 

1.39989 

0.57325 

1.3 

1 

k 

1.400+0.03  x(/-  1) 

2.58649 

0.75869 

1.44236 

1.34749 

1.2 

1 

J 

1.400  +  0.03  v  (/—  1) 

2.72798 

0.95200 

1.50447 

5.56615 

1.1 

1 

1.400  +  0.03  x  (/—  1) 

2.90S96 

1.26615 

1.53330 

131.64487 

19 


Table  8.  BEST  SIMULATION  RESULTS  FOR  MODERATE  SAMPLE  SIZE 


(N  =  50) 


Shape 
param¬ 
eter  a 

{«„}  functional  form 

MSE  of 
actual 
average 
cost 

Esti¬ 

mated 

optimal 

replace¬ 

ment 

time  <£* 

MSE  of 

<j>* 

2.0 

1 

2.225  +  0.06  x  (/—  1) 

1.93167 

0.08290 

1.26858 

0.17087 

1.9 

1 

2.525  +  0.03  x(/-  1) 

1.99418 

0.08626 

1.30783 

0.197S6 

1.8 

1 

(9.00+0.03  x  (/—  l))os 

2.05713 

0.08936 

1.33833 

0.22572 

1.7 

1 

2.225  +  0.03  x(/-  1) 

2.13354 

0.09605 

1.38420 

0.25644 

1.6 

1 

(8.00+0.03  x  (/—  l))0-5 

2.20758 

0.10296 

1.45589 

0.34657 

1.5 

1 

1.400  + 0.03  x(/-  1) 

2.30572 

0.10254 

1.59010 

0.41487 

1.4 

1 

1.400  +  0.03  x  (/—  1) 

2.38234 

0.11150 

1.74888 

0.55554 

1.3 

1 

1.400+0.03  x  (/—  1) 

2.47548 

0.12823 

1.91380 

0.83251 

1.2 

1 

1.400  +  0.03  x  (/—  1) 

2.56903 

0.15299 

2.16169 

3.34231 

1.1 

1 

1.400+0.03  x(/-  1) 

2.68135 

0.20829 

2.35770 

113.78398 
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Table  9. 

Shape 
param¬ 
eter  « 

2.0 

1.9 

1.8 

1.7 

1.6 

1.5 

1.4 

1.3 

1.2 

1.1 


BEST  SIMULATION  RESULTS  FOR  LARGE  SAMPLE  SIZE 


(N  =  250) 


{«„}  functional  form 

Actual 

average 

cost 

Ry 

MSE  of 
actual 
average 
cost 

Esti¬ 

mated 

optimal 

replace¬ 

ment 

time  4 i* 

MSE  of 

i* 

1 

2.525  +  0.06  x  (/—  1) 

1.86401 

0.01816 

1.17291 

0.04966 

1 

2.525  +  0.06  x  (/  -  1) 

1.93160 

0.01  S3 1 

1.20933 

0.05909 

1 

2.525  + 0.03  x(i-  1) 

2.01019 

0.01897 

1.26157 

0.0S345 

1 

2.525  +  0.03  x  (/-  1) 

2.08-192 

0.01918 

1.30963 

0.08491 

1 

(10.0+0.09  x  (/—  l))0-5 

2.17148 

0.01916 

1.41962 

0.14584 

1 

1.725  +  0.03  x  (/—  1) 

2.25858 

0.01997 

1.51001 

0.20026 

l 

1.725  +  0.01  x  (/  —  1) 

2.35111 

0.02039 

1.77728 

0.37160 

1 

1.725  +  0.01  x  (/-  1) 

2.43630 

0.02189 

2.09759 

2.63463 

1 

1.725  +  0.01  x  (/—  1) 

2.51301 

0.02519 

2.60572 

2.22276 

1 

1.725  +  0.01  x(/-  1) 

2.58071 

0.03489 

3.26706 

95.67296 

SMALL  SAMPLE  SIZE  N=10 

1.  a=l .1,  1.2 .  1.7 

V.  2.  a=  1.8 

. .  3.  a »1.9,  2*0 


M0DERATE  SAMPLE  SIZE  N=50 
1.  o-l.l,  1.2, 
N  2.  a=i.6 

3.  a=1.7 

)  ^  4.  a-1.8 

.!  5.  a=1.9 

'  ,  a *=2.0 


NUMBER  OF  REPLACEMENTS  i 


NUMBER  OF  REPLACEMENTS  i 


LARGE  SAMPLE  SIZE  N=250 


i 


1.  a-1.1,  1.2,  1.3,  1.4 

2.  a-1.5 

3.  a  =*1.6 

4.  a-1.7,  1.8 

5.  a=1.9,  2.0 
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IV.  ANALYSIS  OF  SIMULATION  RESULTS 


It  is  clear  from  Tables  7,  8  and  9  that  {<?„}  sequences  with  the  smallest  MSE  when 
the  shape  parameter  a  equals  2.0,  perform  significantly  worse  when  a  is  equal  to  1.1. 
It  is  also  clear  that  for  a  given  a  ,  the  {«„}  sequence  performing  well  for  small  sample 
sizes  will  not  be  the  best  {a,,}  for  large  sample  sizes.  In  this  chapter  we  analyze  in  detail 
how  the  best  performing  {u„}  sequences  differ  for  different  sample  sizes,  values  of  a  for 
the  underlying  system  life  distribution,  and  for  different  costs  C„  C2.  To  make  these 
detailed  comparisons  for  /  =  1,2, ....,  N,  we  will  examine: 

•  c,  ,  the  /til  scheduled  replacement  times  (censoring  ages)  q,  averaged  over  1000  re¬ 
petitions, 

•  P(J,  <  <;,) ,  the  proportion  of  unscheduled  replacements  from  1000  repetitions  at  the 
/th  replacement,  and 

•  R,  ,  the  average  actual  cost  per  unit  time  for  the  first  /  replacements  averaged  over 
1000  repetitions. 

We  begin  by  studying  simulation  results  for  small,  moderate  and  large  sample  sizes  sep¬ 
arately. 

A.  SMALL  SAMPLE  SIZE  N  =  10 

From  Table  7  on  page  19,  there  are  only  three  different  best  performing  {<?„}  se¬ 
quences  for  N  =  10.  We  plot  the  performance  of  all  three  { an }  sequences  for  a  =  2.0, 
1.8,  1.5  and  1.3  in  Figures  5,  6,  7  and  8  respectively.  Figures  5(a)~8(a)  are  the  same; 
they  plot  a,  versus  i,  to  give  a  visual  representation  of  the  three  {a„}  sequences.  Figures 
5(b)~8(b)  plot  q ,  versus  i,  and  the  dotted  line  is  the  optimal  age  of  replacement  <£*. 
Figures  5(c)~8(c)  plot  P(A)<^)  versus  /,  and  the  dotted  line  is  P(A;  <<£*),  the  proba¬ 
bility  of  unscheduled  replacement  under  a  policy  with  optimal  age  of  replacement  0*. 
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Figures  5(d)~8(d)  plot  R,  versus  /  and  the  dotted  line  is  R(0*).  Recall  that  the  first 

scheduled  replacement  time  c,  =  1.0  and  cr,  =  0.0,  so  that  for  each  simulation 

Ci  =  1.0.  In  Figures  5,  6,  7  and  8: 

•  the  curves  marked  “1”  correspond  to  the  {a„}  sequence  of  functional  form 

~1.400  +  a03  x"(7-  1)  which  Performs  best  for  «  =  1.1,  1-2, .  1.7, 

•  the  curves  marked  “2”  correspond  to  the  { an }  sequence  of  functional  form 

_  —  1- — - — — -  which  performs  best  for  a  =  1.8,  and 

(1.72O  +  0.01  x  (/- 1))  ^ 

•  the  curves  marked  “3"  correspond  to  the  {«„}  sequence  of  functional  form 


1.725  +  0.06  x  (/ —  1) 


which  performs  best  for  a  =  1.9,  2.0. 


1.  Shape  Parameter  a  Equals  2.0 

The  average  censoring  time  c,  for  all  three  {a,}  sequences  overestimates  the  op¬ 
timal  replacement  time  0*  =  1.1792,  for  i  =  2,  3,  .....  10  (Figure  5(b)).  Because 
£<  =  4>*i +  a, ,  when  a,  is  large,  the  scheduled  replacement  time  for  the  /th  replacement 
^  is  also  large.  Thus  we  would  expect  that  the  {an}  sequence  with  larger  values  to  \ield 
larger  which  in  turn  would  mean  that  a  higher  proportion  of  replacements  would  be 
unscheduled.  Figures  5(a),  5(b)  and  5(c)  are  consistent  with  this  expectation. 

In  Figure  5(d)  R(0*)  =  1.8024.  This  figure  shows  that  for  all  three  sequences 
the  R,  the  average  cost  per  unit  time  up  to  the  /th  replacement  decreases  with  /,  the 
number  of  replacements.  This  result  is  promising  because  the  goal  of  the  sequential  es¬ 
timation  procedure  is  to  decrease  costs  while  sampling.  Even  though  as  «  ->  oo,  Rn  will 


approach  the  optimal  replacement  cost  R(0*)  with  probability  1.0  [Ref.  8j,  there  is  no 
guarantee  that  Rn  will  decrease  for  the  first  few  observations. 
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Figure  5.  Plot  of  {«„}  performance  when  N  =  10  and  a  =  2.0 


The  ordering  of  the  curves  1,  2  and  3  in  Figure  5(d)  supports  the  intuition  that 
the  {«„}  sequence  which  yields  f,  's  closest  to  0*  will  have  the  smallest  average  costs  per 
unit  time  Rt.  Thus,  curve  3  which  is  closest  to  0*  in  Figure  5(b)  is  also  closest  to 
R(4>*)  in  Figure  5(d). 
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2.  Shape  Parameter  a  Equals  1.8 

Mere  =  1.2032  and  =  1.9470.  Again  all  three  {«„}  sequences  yield 

f,  >  $*,  i  =  2,  3, ....,  10  (see  Figure  6(b)).  However,  comparing  curve  3  in  Figme  6(b) 
with  curve  3  in  Figure  5(b),  we  see  that  the  third  {«„}  sequence  yields  c,  values  closer  to 
4>*:  for  a  =  1.8  than  for  a  =  2.0.  Although  the  third  {<?„}  sequence  clearly  yields  q, 
values  closest  to  0*for  a  =  1.8  (curve  3,  Figure  6(b)),  both  the  second  and  third  {«„} 
sequence  have  nearly  identical  average  actual  cost  per  unit  time  (curves  2  and  3,  Figure 
6(d)).  In  fact,  for  a  =  1.8  the  second  {«„}  sequence  performs  better  than  the  third  {<7„} 
sequence  in  that  the  MSE  for  the  second  {a„}  sequence  is  smaller. 

To  explain  this,  recall  that  at  each  replacement,  we  can  only  observe  the  first  i 
replacements.  During  early  stages  of  sampling  (i.e.,  for  small  /),  the  information  about 
the  system's  life  distribution  will  not  be  enough  to  get  a  very  good  estimates  of  0*. 
Thus  the  estimator  </>*  and  the  scheduled  censoring  time  ^  will  have  large  variances  in 
early  stages  of  sampling.  Also  notice  that  in  Figure  3  on  page  15,  R(t)  to  the  left-hand 
side  of  the  global  minimum  is  steeper  than  on  the  right-hand  side  of  $*.  Thus  slight 
underestimates  of  <f>*  can  increase  costs  much  more  than  corresponding  overestimates 
of  if)*.  This  explains  the  fact  that  when  >  0*  the  {a„}  sequences  with  1,  is  closer  to 
</>*,  may  not  do  as  well  as  {a„}  sequences  with  slightly  large  q,. 

We  conclude  that  when  the  {<?„}  sequence  with  average  censoring  times  q,'s 
closest  to  0*  has  £,></>*,  and  the  differences  the  q,’s  and  <f>*  are  small,  then  this  se¬ 
quence  might  not  yield  the  lowest  average  cost,  and  we  need  to  check  the  variation  in 
the  's  carefully  before  jumping  to  conclusions. 
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SMALL  SAMPLE  SIZE  N  =  10 
(  a  )  (  a  =  1.8  )  (  b  ) 


(c)  (d) 


Figure  6.  Plot  of  {«„}  performance  when  N  -  10  and  «  =  1.8 

3.  Shape  Parameter  a  Equals  1.5 

The  {a„}  sequences  in  Figure  7  perform  similarly  to  those  in  Figure  6.  We  see 
that  the  >  ft*  for  /  =  2,  3,  10  for  all  three  sequences,  and  that  the  three  curves 

of  the  average  censoring  time  are  even  closer  to  ft*  than  when  oc  =  1.8  and  oc  =  2.0 
(Figures  6(b)  and  5(b)).  The  same  line  of  reasoning  used  to  explain  Figure  6  when 
a  =  1.8,  can  be  used  here  to  explain  why  the  average  actual  costs  from  the  three  {c7„} 
sequences  are  virtually  identical  (Figure  7(d))  even  though  the  q,  's  for  the  three  se- 
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quences  are  clearly  separated  (Figure  7(b)).  In  fact,  the  first  {on}  sequence,  with  the 
largest  c,  values  (curve  1,  Figure  7(b))  has  the  smallest  MSE. 


SMALL  SAMPLE  SIZE  N  -  10 
(  o  )  (  a  -  1 .5  )  (  b  ) 


Figure  7.  Plot  of  {«,,}  performance  when  N=  10  and  «  =  1.5 


4.  Shape  Parameter  a  Equals  1.3 

Figure  8(b)  shows  that  for  all  three  {«„}  sequences  g,  <  <p*,  i  —  1,  2, . 10  and 

that  the.first  {«,,}  sequence  (curve  1)  yields  f,  values  closest  to  </>*.  Again  because  R{t ) 
to  the  left  of<£*  is  steeper  than  R{t)  to  the  right  (Figure  3,  curve  2),  when  several  {an} 
sequence  yield  f,  <  <£*  for  /  =  1,2, 10,  then  the  {<?„}  sequence  with  <?/s  closest  to  <£* 
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should  yield  the  lowest  average  actual  costs  per  unit  time  R, .  Thus  it  is  not  surprising 


that  the  first  {«„}  sequence  has  the  smallest  Rt,  i  =  1,  2, 10  and  the  smallest  MSB. 


SMALL  SAMFLE  SIZE  N  = 
(  )  (  a.  =  1.3  ) 


10 


(c) 


(b) 


1  NUMBER  OF  REPLACEMENTS  i 


Figure  8.  Plot  of  {«„}  performance  when  N  =  10  and  *  =  1.3 


B.  MODERATE  SAMPLE  SIZE  N  =  50 

We  use  the  same  approach  to  analyze  the  performance  of  { an }  sequence  with  mod¬ 
erate  sample  size  as  we  used  for  small  sample  size.  In  figures  9,  10  and  11: 

•  the  curves  marked  “1”  correspond  to  the  {«„}  sequence  with  functional  form 

1.400  +  0.03  x  (/-  1)'  Whfch  Performs  b««  «>r  »  -  U  1.2 .  1.5,  . 

•  the  curves  marked  “2”  correspond  to  the  {a„}  sequence  with  functional  form 

(s.oo+o.oiU-or  "'hich  performs  best  for  a=l6’ and 
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•  the  curves  marked  “3"  correspond  to  the  {«„}  sequence  with  functional  form 

—  L — r. — rr  which  performs  best  for  a  =  2.0. 

2.22a +  0.06  x  (/  —  1)  v 

1.  Shape  Parameter  a  Equals  2.0 

Again  in  Figures  9(a),  9(b)  and  9(c)  we  see  the  same  trends  that  were  evident  for 
N  =  10  when  a  =  2.0  (Figures  5(a),  5(b)  and  5(c)  on  page  25).  For  instance,  curve  1 
in  Figure  9(a)  has  the  largest  a,  values,  /  =  1,2, .... ,  50,  in  Figure  9(b)  it  has  the  largest 
c,  values  and  in  Figure  9(c)  it  has  the  largest  proportion  of  unscheduled  replacements. 
Curve  3  has  the  lowest  average  cost  and  it  has  the  smallest  MSE.  Thus  the  {«„}  se¬ 
quence  of  functional  form  ■•  ••—,  — - — —  ,  has  the  best  performance  for  .V  =  50 

1  2.22o +  0.06  x  (i  -  1)  1 


Figure  9.  Plot  of  {«„}  performance  when  N  =  50  and  a  =  2.0 
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2.  Shape  Parameter*  Equals  1.6 

From  Figure  10(d),  we  see  the  same  phenomena  that  occurs  for  N  =  10  with 
a  =  1.8  (Figure  6  on  page  27),  i.e.,  that  although  the  average  scheduled  replacement 
times  c,  for  curve  3  are  closer  to  0*  than  for  curve  2,  curve  2  has  the  lower  average 
cost.  Again  this  can  be  explained  by  the  fact  that  when  c,  >  <f)*  the  variance  of  the 
c,'s  has  a  more  dramatic  effect  on  costs  when  c,  is  close  to  4>*  than  when  c,  is  further 
away.  It  is  interesting  to  note  that  this  effect  of  the  variation  of  $,'s  on  cost  shows  up 
to  same  extent  in  Figure  9  when  a  =  2.0. 


MODERATE  SAMPLE  SIZE  N  =  50 
(a)  (  a  =  1 .6  )  (b) 


Figure  10.  Plot  of  {«„}  performance  when  N  =  50  and  a  =  1.6 
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3.  Shape  Parameter  a  Equals  1.3 

All  three  {u„}  sequences  yield  c,  <  <j>*  for  /  =  1,  2, . 50  (Figure  11(b))  con¬ 

sistent  with  the  observations  made  when  N  =  10  and  a  =  1.3  (Figure  8),  the  first  {«„} 
sequence  (curve  1  in  Figure  11(b))  which  yields  c,  values  closest  to  <f>*,  also  has  the 
smallest  R,'s  (curve  1  in  Figure  11(d))  and  has  smallest  MSE. 


MODERATE  SAMPLE  SIZE  N  =  50 
(  a  )  (  a  =  1.3  )  (  b  ) 


Figure  1 1.  Plot  of  (o„)  performance  when  N  =  50  and  a  —  1.3 
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C.  LARGE  SAMPLE  SIZE  N  -  250 


Figures  12,  13  and  14  show  the  performance  of  {«„}  sequences  as  follows: 

•  the  curves  marked  “1”  correspond  to  the  {«„}  sequence  of  functional  form 

— - — —  which  performs  best  for  a  =  1.3, 

1.723  +  0.01  x  (/  -  1)  r 

•  the  curves  marked  "2"  correspond  to  the  {a,}  sequence  of  functional  form 

(io.oo+b.09  x'(';— which  pcrforms  bcst  for «  “  *-6' and 

•  the  curves  marked  “3"  correspond  to  the  {u„}  sequence  of  functional  form 

— — — 77*  which  performs  best  for  a.  =  2.0. 

2.o2o  +  0.06  x  (<- 1)  * 

1.  Shape  Parameter  a  Equals  2.0 


As  for  the  other  sample  sizes  when  a  =  2.0,  al!  three  {«„}  sequences  yield 


<?,  >  for  /=  2,  3 . .  250.  In  early  stages  of  sampling,  the  average  costs  for  the  second 


{a„}  sequence  (curve  2,  Figure  12(d))  are  lower  than  the  average  costs  for  the  third  {«„} 
sequence  (curve  3,  figure  12(d)).  But  after  about  60  replacements  the  situation  reverses 
and  curve  3  falls  be' nv  curve  ?.  in  Figure  12fd).  From  Figure  12(b),  we  observe  that 
the  average  censoring  times  of  curve  3  are  large  when  /  is  small.  After  about  30  re¬ 
placements,  the  average  censoring  times  of  curve  3  are  the  closest  to  $*.  Because  the 
variance  of  the  estimators  of  0*  should  be  smaller,  after  moderate  number  of  teplace- 
ments,  of  curve  3  (Figure  12(d))  decreases  faster  than  the  other  two  curves  and  fin¬ 
ishes  lower  than  curve  2  (Figure  12(d))  in  the  end. 
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URGE  SAMPLE  SIZE  N  =  250 
(  a  )  (  a  =  2.0  )  (  b  ) 


Figure  12.  Plot  of  {«„}  performance  when  N  =  250  and  a  =  2.0 


2.  Shape  Parameter  a  Equals  1.6 

Although  the  second  and  third  {a„}  sequences  are  fairly  different  (curves  2,  3  in 
Figure  13(a)),  their  average  costs  are  about  the  same.  Before  these  two  {«„}  sequences 
intersect,  the  proportion  of  uncensored  observations  (i.e.,  unscheduled  replacements)  is 
higher  for  the  second  {«„}  sequence  than  the  third  (curves  2  and  3,  Figure  13(c)).  After 
intersection  the  second  {<?„}  sequence  has  a  higher  proportion  of  unscheduled  replace¬ 
ment.  Thus  after  a  moderate  number  of  replacements,  it  would  seem  that  about  the 
same  amount  of  information  about  the  system's  life  distribution  and  is  gathered  using 
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either  {an}  sequence,  so  that  we  could  expect  the  average  costs  to  be  about  the  same. 
However,  after  a  large  number  of  replacements,  the  variance  of  censoring  times  become 
small  and  the  {a„}  sequence  yielding  c,  's  the  closest  to  the  optimal  replacement  time 
(curve  3)  has  the  lowest  average  costs. 

Although  the  third  { a„ }  sequence  (curve  3  in  figure  12(d))  gives  the  lowest. 
our  simulation  results  reveal  that  the  variance  of  actual  costs  per  unit  time  are  larger 
than  for  the  second  {a*}  sequence.  Our  policy  is  to  choose  the  sequence  with  the 
smallest  .VISE  as  the  best  performer.  Thus  we  conclude  that  the  second  {<7„}  sequence 
(curve  2),  is  the  best  performer. 


LARGE  SAMPLE  SIZE  N  =  250 
(  a  )  (  a  =  1.6  )  (  b  ) 


CO  (d) 


Figure  13.  Plot  of  {«„}  performance  when  N  =  250  and  «  =  1.6 


3.  Shape  Parameter  a  Equals  1.3 

Since  the  second  and  third  {«„}  sequences  (curves  2  and  3  in  Figure  14(b))  yield 
<*;,  <  (/>*  for  i  -  1,2,  250,  we  expected  the  sequences  with  closer  to  <j)*  to  have 

the  lowest  cost.  This  can  be  seen  with  curves  2  and  3,  we  see  in  Figure  14(b).  The  first 
sequence  (Curve  1  in  Figure  14(b))  yields  f/s  closest  to  <£*.  Thus  it  is  not  surprising 
to  find  out  that  the  first  sequence  (curve  1  in  Figure  14(d))  has  the  lowest  average  cost 
and  the  smallest  MSE. 


LARGE  SAMPLE  SIZE  N  =  250 
(  a  )  (  a  =  1.3  )  C  b  ) 


Figure  i4.  Plot  of  {«„}  performance  when  N  =  250  and  a  =  1.3 
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D.  SAMPLE  SIZE  CHANGES 


In  the  last  section  for  each  fixed  sample  size  small,  moderate  and  large,  we  notice 
that  as  the  shape  parameter  a  decreases  from  2.0  to  1.1,  the  best  performing  {u„}  se¬ 
quences  tend  to  have  larger  values.  This  means  that  when  the  system's  life  distribution 
is  close  to  the  exponential  distribution,  the  seq.  .ntial  estimation  piocedure  needs  to  use 
{«„}  sequences  with  larger  values  to  get  smaller  actual  costs  per  unit  time.  Now  we  ex¬ 
amine  how  the  {a,,}  sequences  perform  as  sample  size  \  aries  for  each  fix  u.  We  picked 
some  of  the  best  performing  { an }  sequences  from  Tables  7,  8  and  9,  and  plotted  these  {u„} 
sequences  in  Figure  15.  The  short  curves  are  the  best  performing  {<?„}  sequences  for 
small  sample  size  .V  =  10,  the  moderate  curves  are  the  best  performing  {<?„}  sequences 
for  moderate  sample  size  A'  =  50  and  the  long  curves  are  the  best  performing  {«,}  se¬ 
quences  for  large  sample  size  .V  =  250. 
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In  the  bottom  plot  of  Figure  15,  the  {u„}  sequence  were  the  same  for  a  =  1.4,  1.3, 
1.2  and  1.1.  We  already  know  that  as  the  system's  life  distribution  approaches  the  ex¬ 
ponential  distribution  (i.e.,  as  the  shape  parameter  a.  decreases  to  1.0)  the  sequential  es¬ 
timation  procedure  performs  best  with  {«-„}  sequences  with  larger  values.  Thus  the 
condition  that  four  different  a  values  have  the  same  best  performing  sequences,  is  due 
to  the  fact  that  {<?„}  sequences  with  large  enough  values  were  not  tried  in  the  simulation. 

From  Figure  15,  for  each  a,  the  best  {a„}  sequences  for  N  =  10  are  larger  than  the 
best  sequences  for  N  =  250.  However  best  {«,,}  sequences  for  N  =  50  sometimes 
intersect  with  the  best  {a„}  sequences  for  A’  =  250.  Generally  speaking,  the  relationship 
between  three  best  {<?„}  sequences  is  not  consistent  for  different  a  values.  In  early  stages 
of  sampling  before  much  data  been  collected,  the  variation  of  censoring  ages  is  very 
large.  From  Figure  3  on  page  15,  if  the  censoring  time  is  less  than  the  optimal  re¬ 
placement  time  <f>*,  then  the  censoring  time  will  be  on  the  steeper  side  of  the  curve.  This 
means  the  censoring  time  will  have  a  large  expected  average  cost.  Thus  we  would  rather 
have  an  { an }  sequence  with  large  values  so  that  the  qt  's  tend  to  be  larger  than  0*  to 
reduce  the  risk  of  higher  costs  when  <  0*. 

When  /  becomes  large,  we  have  many  observations  with  more  information  about  the 
system's  life  distribution  which  can  be  used  to  better  estimate  0*,  thereby  reducing  the 
variance  of  the  c/s.  Thus  the  {«„}  sequences  giving  average  censoring  times  q,  closest  to 
the  optimal  replacement  time  0*  performed  best.  In  general  {a„}  needs  to  decrease  as 
n  increases.  However  for  distributions  that  are  closer  to  exponential,  0*  is  larger  and 
the  estimators  0*  tend  to  be  less  than  0*  ,  we  need  { an }  sequences  with  large  values  that 
decrease  slowly.  For  large  a  values,  0*  tends  to  overestimate  the  optimal  replacement 
time,  thus  the  {«„}  sequence  which  drop  to  zero  faster  or  have  light  tails,  achieves  lower 
cost. 
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E.  COST  RATIO  CHANGES 


1.  Optimal  Replacement  Time 

We  recall  that  the  expected  long  run  costs  per  unit  time  R(t)  is  dependent  on  the 
unscheduled  and  scheduled  maintenance  cost,  the  system's  life  distribution  and  the 
scheduled  censoring  time  /: 


R(t)  = 


C,  x  F(i)  +  C2  x  S(t) 


I  S(x)dx 
Jn 


where  C,  and  C2  are  unscheduled  and  scheduled  replacement  cost  respectively.  By  ex¬ 
pressing  R{t)  as 


C, 


+S(  7) 

R(t)  =  C2  x  — ^ - 

'  !  S(x)  d  x 


C, 


we  see  that  the  optimal  replacement  time  0*  is  a  function  of  cost  ratio  .  Thus  for 
a  given  system  life  distribution,  as  long  as  the  ratio  is  constant,  then  </>*  is  the  same  no 
matter  how  large  C,  and  C2  are.  Thus  we  cun  study  the  effect  of  different  costs  on  the 
sequential  estimation  procedure  by  changing  the  cost  ratio.  Table  10  shows  optimal 
replacement  times  for  several  cost  ratios.  These  were  found  using  the  same  method  used 
to  construct  Table  5  on  page  16. 
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Table  10.  COMPARISON  OF  OPTIMAL  REPLACEMENT  TIMES  BY  DIF¬ 
FERENT  COST  RATIO  WITH  E(X;)=2.0 


Shape 

Parameter  a 

- ! 

Optimal  Replacement  Time  <£* 

C,  =  2.0 

C2  =  1.0 

B 

C,  =  8.0 

C2  =  1.0 

C,  =  1 0.0 

C2  =  1.0 

2.0 

2.4638243 

1.1792297 

0.S876661 

0.737S923 

1.9 

2.7862625 

1.1846266 

0.9107760 

0.6948394 

1.8 

2.9159756 

1.2031918 

0.S641901 

0.7974927 

1.7 

3.3440571 

1.2661085 

0.9595973 

0.8142039 

1.6 

3.6384935 

1.3114500 

0.9868178 

0.7605762 

1.5 

4.8073025 

1.4643221 

1.0388746 

0.9282054 

1.4 

7.6937S76 

1.6853428 

1.0966682 

0.9166338 

1.3 

10.2272406 

2.3491287 

1.2848225 

1.0260086 

1.2 

18.5055695 

3.7472544 

1.7184162 

1.4291582 

1.1 

21.9620056 

12.9815130 

4.2228003 

3.S461828 

From  Table  10  for  each  fixed  ratio  costs,  the  optimal  replacement  time  </>*  in¬ 
creases  as  the  shape  parameter  a  decreases  from  2.0.2  From  Table  10  we  also  see  that 

C 

for  fixed  a  as  the  cost  ratio  ~  increases,  the  optimal  replacement  time  <j>*  decreases. 

c2 

This  means  that  when  the  unscheduled  replacement  cost  is  too  large,  we  don't  want  to 
risk  unscheduled  replacement. 


2  Note  that  in  Table  10,  when  C,  =  8.0,  the  optimal  replacement  time  cf>*  value  for  a  = 
1.8  is  less  than  the  <p*  value  for  a  =  1.9.  The  same  is  true  when  C,  =  10.0,  a  =  1.6  and  1.4. 
This  is  caused  by  the  specific  series  of  Pseudo  Random  Numbers  (PRN)  used  to  simulate  <p*. 
When  the  seed  for  the  PRN  generator  changes,  the  simulation  results  for  these  situations  should 
change  also. 
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2.  Effect  on  Performance 

Let  a  =  1.8,  the  scheduled  replacement  cost  C2  =  1.0,  and  let  the  unscheduled 
replacement  cost  C,  have  different  values  2.0,  5.0,  8.0  and  10.0.  To  examine  the  per- 

(7 

formance  of  the  {an}  sequence  for  different  ratios,  we  simulated  the  sequential  esti¬ 
va 

mation  procedure  with  each  of  the  { a „}  sequences  from  Table  6  on  page  18.  Again,  the 
best  performing  {crn}  sequences  were  chosen.  They  are  summarized  in  Table  11. 


Table  11.  BEST  PERFORMING  {«„}  FOR  SHAPE  PARAMETER  a  =1.8 


Unscheduled 

and 

Scheduled 

Replacement 

Cost 

{«„}  Functional  Form 

_ 

A^=  10 

A  «  50 

N  =  250 

C,  =  2.0 

C:=  1.0 

1 

1 

1 

1.400  +  0.03  x(/-  1) 

1.400  +  0.03  x  (/—  1) 

1.725  +  0.01  x  (/—  1) 

C,  =  5.0 

C=  1.0 

1 

1 

1 

1.725  +  0.01  x  (/—  1) 

(9.0+0.03  x  (/—  1))0  J 

2.525  +  0.03  x  (/  —  1) 

C,  =  8.0 

C,=  1.0 

1 

1 

1 

2.525  +  0.03  x(i-  1) 

2  525  +  0.15  x  (/  -  1) 

2.525  +  0.15  x  (/-  1) 

C,=  10.0 
C=  1.0 

1 

1 

1 

2.525  +  0.03  x  (/  —  1) 

2.525  +  0.15  x  (/  -  1) 

2.525  +  0.15  x  (/  -  1) 

From  Table  11,  we  observe  that  for  C,  =  8  and  10,  the  same  {a,,}  sequences 

perform  best  for  sample  sizes  N  =  10,  50  and  250.  This  could  be  because  the  cost  ra- 
C 

tios  ~  =  8  and  10  are  close,  or  it  might  be  caused  by  the  fact  that  we  did  not  test 
L2 

enough  {aj  sequences. 
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Figure  16.  The  {«„}  sequence  performances  for  a  =  1.8 
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F.  COMPARISON  WITH  FREES  &  RUPPERT'S  PROCEDURE 


There  have  been  two  previous  attempts  to  construct  nonparametric  sequential  esti¬ 
mators  of  0*:  Bather  [Ref.  11:  pp.  253-266]  and  Frees  &  Ruppert  [Ref.  12: 
pp.  650-661].  Bather's  procedure  was  not  fully  sequential  in  that  periodically,  the  sys¬ 
tem  was  allowed  to  operate  until  failure  regardless  of  how  long  it  took.  Frees  and 
Ruppert  use  stochastic  approximation  techniques  to  construct  a  fully  sequential  esti¬ 
mator  of  <f>*  (i.e.,  c,  <  oo  for  /=  1, 2, ....).  In  this  section  we  compare  the  performance 
of  the  sequential  procedure  described  in  Chapter  3  with  the  perform.mce  of  Frees  and 
Ruppert's  procedure.  Because  Frees  and  Ruppert  do  not  estimate  directly,  we  can 
not  compare  the  procedures  using  mean  square  error  (MSE)  of  the  estimator  o f  0*. 
We  can  however  compare  MSE,  the  squared  difference  of  actual  costs  per  unit  time  from 
/?(</>*)  averaged  over  the  repetitions  of  each  simulation,  for  both  procedures.  Because 
the  ultimate  goal  is  to  reduce  costs,  this  is,  in  fact,  the  more  meaningful  comparison. 

In  their  Monte-Carlo  study,  Frees  and  Ruppert  took  the  unscheduled  replacement 
cost  C,  =  5.0  and  the  scheduled  replacement  cost  C2  =  1.0.  The  system's  life  distrib¬ 
ution  is  Weibull  with  shape  parameter  a  =  2.2  and  scale  parameter  X  =  0.5.  This  dis¬ 
tribution  has  expected  value  1.7712  and  standard  deviation  0.8499.  The  optimal 
replacement  time  4>*  =  0.99505  and  optimal  cost  =  1.904.  [Ref.  12:  p.  656] 

Frees  and  Ruppert  update  their  estimator  of  </>*  after  every  two  observations.  They 
chose  to  look  at  the  performance  of  their  estimator  after  10,  50  and  250  updates  to  re¬ 
flect  small,  moderate  and  large  sample  sizes  respectively.  Thus  after  10,  50  and  250  up¬ 
dates,  they  used  20,  100  and  500  observations  respectively.  Table  12  on  page  44  shows 
the  performance  of  cost  estimator  from  Frees  and  Ruppert's  sequential  estimation  pro¬ 
cedure. 
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In  order  to  compare  methods,  we  use  the  same  simulation  setting  as  Frees  and 
Ruppert's,  C,  =  5.0,  C2  =  1.0,  the  system's  life  distribution  is  Weibull  with  shape  pa¬ 
rameter  a  =  2.2  and  scale  parameter  A  =  0.5,  and  we  select  a  { an }  sequence  of  the 
functional  form  .  This  { a„ }  sequence  roughly  corresponds  to  the  sequence 

used  by  Frees  and  Ruppert.  The  performance  of  our  cost  estimator  is  shown  in  the 
Table  13. 


Table  12.  THE  PERFORMANCE  OF  COST  ESTIMATOR  FROM  FREES  AND 
RUPPERT'S  PROCEDURE 


Number  of  updates 

10 

50 

250 

Number  of  observations 

20 

100 

500 

Average  cost 

2.268 

2.159 

2.053 

Variance  of  average  cost 

0.1851 

0.0408 

0.0096 

MSE  of  average  cost 

0.3176 

0.1058 

0.031S 

Table  13.  THE  PERFORMANCE  OF  COST  ESTIMATOR  FROM  OUR  SE¬ 


QUENTIAL  ESTIMATION  PROCEDURE 


Number  of  updates 

10 

50 

250 

Number  of  observations 

10 

50 

250 

Average  cost 

2.168 

2.005 

1.952 

Variance  of  average  cost 

0.6325 

0.0971 

0.0158 

MSE  of  average  cost 

0.7022 

0.1073 

0.0181 

Number  of  updates 

20 

100 

500 

Number  of  observations 

20 

100 

500 

Average  cost 

2.086 

1.983 

1.940 

Variance  of  average  cost 

0.2874 

0.0415 

0.0087 

MSE  of  average  cost 

0.3205 

0.0477 

0.0100 
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Recalling  that  Frees  and  Ruppert's  need  two  observations  for  each  update  of  their 
estimator  of  </>  and  we  only  need  one  observation  to  update  our  estimator  of  <£*.  The 
results  of  our  simulation  after  20,  100  and  500  observations  arc  comparable  to  Frees  and 
Ruppert's  simulation  with  10,  50  and  250  updates  of  </>*  respectively.  Comparing 
Table  13  to  Table  12  with  the  same  number  of  updates,  the  variance  of  our  simulated 
average  costs  are  slightly  higher  than  Frees  and  Ruppert's  simulation  results.  However: 

•  Our  average  costs  are  much  lower  than  Frees  and  Ruppert's  average  costs. 

•  For  small  sample  size  (20  observations),  for  this  {a,}  sequence  ,  the  \1SE  for  both 
procedures  are  about  the  same.  For  large  sample  size  the  MSE  of  our  procedure 
is  smaller  than  that  of  Frees  and  Ruppert,  even  using  half  the  number  of  observa¬ 
tions. 

•  This  choice  of  {<?„}  gives  the  best  performance  for  Frees  and  Ruppert's  procedure 
and  was  chosen  knowing  that  the  underlying  distribution  is  Weibull  (a  =2.2, 
z  =  0.5).  It  is  probable  that  a  different  choice  of  {«?„}  would  yield  better  results  for 
our  sequential  procedure. 

•  The  intermediate  estimators  of  Frees  and  Ruppert's  sequential  estimation  proce¬ 
dure  are  complicated  to  calculate. 

Thus  we  may  conclude  that  our  sequential  estimation  procedure  for  age  replacement 
policies  is  competitive  to  Frees  and  Ruppert's  procedure,  and  is  in  many  ways  better. 
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V.  GRAPHICAL  DETERMINATION  OF  THE  SCHEDULED 
REPLACEMENT  TIME 

A  graphical  method  can  be  used  to  estimate  the  optimal  age  replacement  policy 
based  on  observing  a  sample  of  Lid.  system  lifetimes  [Ref.  3:  pp.  113-115].  This 
method  obtains  the  estimate  of  4>*  from  a  scaled  total  time  on  test  plot.  The  advantage 
of  this  method  is  that  you  can  see  immediately  how  sensitive  the  estimate  of  0*  is  to 
change  in  the  cost  ratio.  This  method  can  be  extended  to  estimate  <£*  based  on  data 
from  our  sequential  estimation  procedure.  In  this  chapter,  we  describe  the  graphic  de¬ 
termination  method  for  i.i.d.  data  and  apply  it  to  the  22  D9G-66A  Caterpillar  Tractor 
Engines  (this  data  was  used  as  an  example  in  Chapter  2).  We  also  indicate  how  this 
method  can  be  extended  to  our  sequential  estimation  procedure. 

A.  SCALED  TOTAL  TIME  ON  TEST  PLOT 

Plotting  the  data  is  often  the  first  step  to  unlocking  information  contained  in  data 
about  the  underlying  model.  For  example,  a  plot  of  the  empirical  distribution  function 
contains  information  about  the  probability  density.  If  the  empirical  distribution  looks 
concave,  then  the  density  may  be  decreasing.  However,  this  plot  is  scale  dependent  so 
that  the  perceived  shape  of  the  plot  depends  on  the  choice  of  plotting  scale. 

A  total  time  on  test  plot  provides  information  about  the  failure  rate.  Analyzing 
failure  data,  it  is  often  the  failure  rate  which  is  of  chief  interest.  For  example,  if  the 
failure  rate  is  constant  or  decreasing,  we  know  that  we  should  not  adopt  a  scheduled 
replacement  policy,  since  an  old  system  in  this  case  is  actually  “better”  than  a  new  sys¬ 
tem.  Consider  an  experiment  where  n  systems  system  lifetimes  A'„  X2,  ....,  Xn  are  ob¬ 
served.  Let  n(u)  be  the  number  of  systems  that  survive  to  age  u,  then  the  total  time  on 
test  to  age  x  is  defined  to  be 
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Let  *„>  ^  <  ....  <  be  the  order  statistics  of  Xxt  A'lt .....  Xn,  then 


r  *» 

T(X{1))  =  n(u)  du 


=  nX(i)  +  (n—  1)(A'(2)  -  X{{))  +  ....  +  (//  -  i  +  l)(A"(i)  - 


is  the  total  time  on  test  to  age  X(i).  We  call 


7'(.v) 


TO 


the  scaled  total  time  on  test  at  age  .r.  A  plot  of  versus  the  empirical  distrib- 

ution  F„(X{0)  =  —  for  /  =  1,2, n  is  called  the  scaled  total  time  on  test  plot.  This  plot 
provides  information  about  the  failure  rate.  If  the  total  time  on  test  plot  is  strongly 
concave,  then  there  is  evidence  that  the  underlying  distribution  is  IFR  and  a  scheduled 
age  replacement  policy  makes  sense. 

B.  OBTAINING  THE  AGE  REPLACEMENT  TIME  FROM  THE  PLOT 

Let  C,  be  the  cost  (in  dollars)  of  an  in-service  tractor  engine  failure  replacement  and 
C2  be  the  cost  (in  dollars)  of  a  scheduled  tractor  engine  replacement,  where  C,  >  C2  . 
After  22  observations,  <f>*  can  be  estimated  from  the  tractor  engine  failure  ages  using  the 

procedure  described  in  Chapter  2  with  a,  -  t,  =  oo,  /  =  1,2 . 22.  It  can  also  be 

estimated  using  the  scaled  total  time  on  te  >.  plot.  Table  14  on  page  48  shows  the  or¬ 
dered  actual  failure  age  of  22  D9G-66A  Caterpillar  Tractor  Engines,  and  the  empirical 

distribution  of  tractor  engine  failures  is  plotted  in  Figure  17  on  page  49. 

—  C 

To  estimate  <£*  plot  the  point  -rz — brr  on  the  x-axis  and  draw  a  tangent  line  to 

(Ct  —  C2) 

the  scaled  total  time  on  test  plot  through  this  point.  Then  drop  down  to  the  x-axis  from 


Table  14.  THE  ORDERED  AGE  OF  FAILURE  TRACTOR  ENGINES 


Order  sequence 
of  failure  age 

Tractor's  number 

Date  of  failure 

Age  in  hours  when 
failure  occurred 

1 

10 

5-8-67 

2690 

2 

20 

7-30-71 

3286 

3 

3 

10-11-71 

3826 

4 

8 

8-8-69 

4394 

5 

21 

1-21-72 

4S15 

6 

2 

4-16-70 

5085 

7 

16 

4-15-65 

5161 

8 

12 

6-28-65 

5278 

9 

19 

4-3-67 

5556 

10 

5 

6-1-70 

6052 

11 

22 

10-31-69 

6150 

12 

11 

3-26-68 

6259 

13 

6 

6-7-71 

6367 

14 

13 

8-1-66 

6378 

15 

14 

9-14-66 

6385 

16 

15 

8-3-66 

6578 

17 

17 

10-26-66 

6717 

18 

18 

11-1-67 

6869 

19 

7 

8-10-70 

7774 

20 

1 

6-16-71 

8230 

21 

9 

9-21-70 

10517 

22 

4 

5-8-72 

10950 

the  point  of  tangency  closest  to  the  value  -jj*  when  i0  =  1,  2, ....,  n.  Then  X(,o)  will  be  the 
estimated  optimal  replacement  age.  A  straight  forward  geometric  argument  shows  that 
estimating  <£*  from  the  scaled  total  time  on  test  plot  is  identical  to  the  estimate  obtained 
by  minimizing  R{T)  [Ref.  3:  p.  114]. 


EMPIRICAL  DISTRIBUTION  Fn 


Figure  17.  The  empirical  distribution  and  the  total  time  on  test  plot  for  Caterpillar 
Tractor  Engine  data 


C.  THE  TOTAL  TIME  ON  TEST  PLOT  OF  TRACTOR  FAILURE  ENGINES 

The  bottom  plot  in  Figure  17  is  the  scaled  total  tir'e  -n  test  plot  when  the  un¬ 
scheduled  replacement  cost  C,  =  200  and  the  scheduled  replacement  cost  C2  =  100. 

-C2 


On  the  horizontal  axis  we  can  find  the  point 


=  -1.0.  We  draw  a  tangent  line 

(C,  —  L2) 


to  the  scaled  total  time  on  test  plot  through  this  point  and  drop  down  to  the  x-axis  from 

I  6 

the  point  of  tangency,  we  find  the  value  —  m  0,273  ^  Then  the  order  statistic  J(6> 
(i.e.,  5085  hours),  is  the  estimated  optimal  replacement  time.  Changing  the  unscheduled 
replacement  cost  C,  to  400  and  700,  and  holding  the  scheduled  replacement  cost  C2  at 
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100,  from  the  scaled  total  time  on  test  plots  in  Figure  18  on  page  50,  we  see  that  the 
estimated  optimal  replacement  age  is  4394  hours  and  2690  hours  respectively.  The  re¬ 
sults  obtained  from  this  figure  are  consistent  with  the  optimal  replacement  ages  esti¬ 
mated  using  the  sequential  estimation  procedures  in  Table  4  on  page  11.  The  value  of 
the  estimator  using  the  scaled  total  time  on  test  plot  depends  upon  how  closely  the  es¬ 
timated  distribution  from  the  sample  approximates  the  true  underlying  life  distribution. 


SCALED  TOTAL  TIME  ON  TEST  PLOT 


Figure  18.  The  total  time  on  test  plot  for  Caterpillar  Tractor  Engine  data 
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We  note  that  this  procedure  can  be  adapted  to  the  sequential  estimation  procedure 
[Ref.  8]  by  replacing  the  empirical  distribution  as  an  estimator  off,  by  the  product-limit 
estimator  of  F  described  in  Chapter  2  and  defining 


T(x) 

AZin)) 


rx 


S(u )  du 


S(u)  du 


where 


S(u) 


=  f]  (  -  ~  — ■  V 

V  /.*--/+  1  J 


T{Z, ) 

is  the  product-limit  estimator  of  S  =  1  —  F.  Then  plot  versus  F(ZW)  -  1  -  S(Z(0) 

*  yz-w 

to  get  the  scaled  “total  time  on  test”  plot.  Following  the  same  steps  outlined  for  the 
-Ca 

u.a.  case,  plot 


c,-c2 


on  the  x-axis  and  draw  the  tangent  line  to  the  scaled  total  time 


on  test  plot  through  the  point 


-C, 


.  The  estimator  of  <p*  is  F*l(.v0)  where  x0  is  the 


c,-c2 

value  on  the  x-axis,  below  the  point  where  the  tangent  line  intersects  the  scaled  total 
time  on  test  plot. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Under  an  age  replacement  policy,  a  system  is  replaced  at  failure  or  after  being  in 
service  for  t  units  of  time,  whichever  comes  first.  The  time  t  is  called  the  scheduled 
censoring  time  or  the  scheduled  replacement  time.  An  is  optimal  replacement  time  0* 
achieves  the  smallest  long  run  expected  cost.  An  important  problem  is  the  estimation 
of  0*  when  the  form  of  the  system's  underlying  life  distribution  is  unknown.  We  show 
by  example  and  through  simulation  that  substantial  cost  savings  can  be  effected  using 
the  sequential  estimation  procedure  for  0*describcd  in  Chapter  2. 

An  important  part  of  our  analysis  is  to  provide  practical  guidelines  for  choosing 
q,  =  0,-i  +  the  schedv  ed  replacement  time  to  be  used  for  the  <th  system,  where  0*,  is 
the  estimator  of  0*  based  on  (/  -  1)  replacements,  and  {a„}  is  a  fixed  sequence  of  con¬ 
stants.  From  our  analysis,  in  the  early  stages  of  sampling  we  would  rather  choose  large 
values  for  a,  so  that  tends  to  be  greater  than  0*.  If  the  scheduled  replacement  times 
q,  are  too  small  in  the  early  stages,  then  most  of  the  system  life  times  are  censored.  Thus, 
very  little  information  is  collected  about  the  underlying  distribution.  This  leads  to  poor 
estimates  for  0*.  After  a  large  number  of  replacements,  the  estimates  of  the  underlying 
distribution  and  therefore  0*  are  relatively  good.  A  good  age  replacement  policy  will 
keep  the  scheduled  censoring  time  closely  approximate  the  optimal  replacement  time 
0*.  Therefore  the  best  performing  {<?„}  sequences  have  smaller  a,  for  large  /. 

The  proper  choice  of  {a„}  sequence  is  a  very  difficult.  The  best  choice  of  an  { a„ } 
sequence. fot;  sequential  estimation  of  age  replacement  policy  depends  on  the  system's 
underlying  life  distribution  which  is  generally  unknown,  the  costs  C„  C2,  and  the  sample 
size  A’  at  which  sampling  stops.  From  oqr  simulation,  we  see  that  when  the  underlying 
distribution  is  close  to  exponential  the  {<?„}  sequences  should  have  larger  values  and 
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should  not  decrease  as  fast  as  those  whose  failure  rate  increases  more  rapidly.  For  large 


sample  sizes  {c7„}  sequences  should  decrease  faster  in  early  stages  than  small  sample  sizes. 


In  addition  as  the  cost  ratio 


£l 

Q 


increases,  the  best  performing  {«„}  sequence  have 


smaller  values. 


Future  research  might  be  concerned  with  the  following  subjects 

•  Sequential  estimation  \  hen  the  underlying  life  distribution  F  comes  from  a 
parametric  family  of  distributions. 

•  How  to  use  the  data  to  estimate  the  {an}  sequence. 

•  If  an  imperfect  repair  model  is  permitted,  then  the  next  failure  depends  on  system 
history,  failure  distributions  are  not  independently  and  identically  distributed 
(/./,</.),  how  to  change  our  Sequential  Estimation  Procedure  to  fit  new  situation. 

•  Minimizing  long  run  expected  costs  per  unit  time  gives  the  same  optimal  <p*f or 
C,  =  2,  Q  =  1  and  costs  for  C,  =  2,000,000,  C2  =  1,000.000.  Clearly  minimizing 
long  run  expected  costs  is  not  appropriate  under  all  circumstances.  Other  cost 
functions  need  to  be  considered. 
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